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Abstract 

We consider a fermion gas on a star graph modeling a quantum wire junction and 
derive the entanglement entropy of one edge with respect to the rest of the junction. 
The gas is free in the bulk of the graph, the interaction being localized in its vertex 
and described by a non-trivial scattering matrix. We discuss all point-like interactions, 
which lead to unitary time evolution of the system. We show that for a finite number of 
particles N, the Renyi entanglement entropies of one edge grow as In N with a calculable 
prcfactor, which depends not only on the central charge, but also on the total transmission 
probability from the considered edge to the rest of the graph. This result is extended to 
the case with an harmonic potential in the bulk. 
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1 Introduction 



Quantum field theory on graphs attracted recently much attention mainly in relation with the 
study of the transport properties of quantum wire networks. Different frameworks [7]-pU] 

have been developed to investigate the phase diagram and the conductance of these structures. 
Despite of the fact that the universal properties in the bulk are described by the well known 
Luttinger liquid theory, the different boundary conditions at the junctions lead to exotic phase 
diagrams [lOl [TTl [121 [151 [161 [T71 [HI [19] whose degree of universality is not completely understood 
and is still under investigation. The results, concerning the charge transport, confirm that the 
conductance properties of the quantum wire networks are strongly affected by the boundary 
conditions as well. 

In the present paper we analyze another physical quantity - the entanglement entropy of one 
edge of the junction with respect to all the others edges. Lots of studies on the entanglement 
properties of many-body systems in the last decade have unveiled new (universal) features of 
these systems and somehow put their global understanding on a deeper level (see e.g. the 
reviews [21] )• In particular, von Neumann and Renyi entanglement entropies of the reduced 
density matrix pA of a subsystem A turned out to be particularly useful for ID systems. Renyi 
entanglement entropies are defined as 

5^°) = ^^InTrp^. (1.1) 
1 — a 

For a —7- 1 this definition gives the most commonly used von Neumann entropy S^^^ = 
—TipA In Pa, while for a — )■ cxd is the logarithm of the largest eigenvalue of pA also known 
as single copy entanglement [22]. Furthermore, the knowledge of the S*^"^ for different a char- 
acterizes the full spectrum of non-zero eigenvalues of pa |23j . 

One of the most remarkable results is the universal behavior displayed by the entanglement 
entropy at ID conformal quantum critical points, determined by the central charge [21] of the 
underlying conformal field theory (CFT) [251 EEl [271 EH]- For a partition of an infinite ID 
system into a finite piece A of length i and the remainder, the Renyi entanglement entropies 
for i much larger than the short-distance cutoff a are 

S^'^^ = ^(l + -]\n- + c^, (1.2) 
\ a / a 

where c is the central charge and Cq, a non- universal constant. 

Given the importance of this result (and also many others not mentioned here) for homo- 
geneous systems, it is natural to wonder whether in the case of junctions the entanglement 
entropies can share some light on the universality and on the relevance of the parameters defin- 
ing the junction. Previous studies in the subject [29l [301 ED [321 EHl [311 ESI |361 137] have been 
limited to the case of only two edges (i.e. an infinite line with a defect) and most often per- 
formed for lattice models. These results provide a strong evidence that the logarithmic behavior 



in Eq. (1.2) remains valid even in the presence of defects, but the prefactor does not depend 
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only on the central charge of the bulk CFT when the defect is a marginal perturbation (in 
renormalization group sense) as it is known to happen for free fermions [1]. 

In order to tackle the problem of entanglement in a junction with an arbitrary number of 
wires, we use the recently developed systematic framework [3HI ES] for calculating the bipartite 
entanglement entropy of spatial subsystems of one- dimensional quantum systems in continuous 
space. We only consider a free fermion gas in bulk in which the junction introduces a marginal 
perturbation. The junction boundary conditions define a specific scattering matrix S, encoding 
all possible point-like interactions in the vertex which give rise to unitary time evolution. 
Focussing on the scale invariant case, we show that for a finite number of particles and 
for edges of equal length L, the Renyi entanglement entropies of any of the edges grow as In A^. 
Oppositely to the case in the absence of the point-like interaction (i.e. Eq. ( 1.2[ )), the prefactor 
of this logarithm does not depend only on the central charge, but also on the total transmission 
probability (1 — |§iip) from the considered edge i to the rest of the graph. Some of the results 
presented here have been anticipated in the short communication [3H]- We show also that the 
presence of an external harmonic potential in the bulk (acting identically on all edges) does not 
alter this result. 

The paper is organized as follows. In the next section we describe the basic features of the 
model and the scattering matrices generated by the point-like interactions at the junction. We 
discuss in detail the scale invariant case and derive the two-point correlation function. The 
entanglement entropy, associated with this system, is analytically computed in section 3. In 
section 4 we extend our considerations, adding a harmonic potential in the bulk. Section 5 is 
devoted to the conclusions and the discussion of some further developments in the subject. 



2 Schrodinger junction 
2.1 The general setting 

In this section we consider a gas of spinless fermions on a quantum wire junction. We 
consider only the ground-state of such system and we refer to it as the ground-state of 
particles, having in mind that the particles are the original fermions and not the excitations 
above the ground-state (that are usually referred as particles in field theory literature). A 
simple model, describing the junction, is represented by a star graph T with M edges of finite 
length L, as shown in Fig. [l] Each point P in the bulk of F is parametrized by {x,i), where 
< X < L is the distance of P from the vertex V of the graph and i labels the edge. We 
assume that in the bulk {x ^ and x ^ L) the gas is free and is described by the Schrodinger 
field ipiit,^), which satisfies 

[^9t + ^dl^Ut,x) = (2.1) 

and standard equal-time canonical anticommutation relations. The only non-trivial interactions 
are localized in the vertex V of F and are encoded in boundary conditions at x = 0. These 
conditions are fixed in turn by imposing that the bulk Hamiltonian —d^ admits a self-adjoint 
extension on the whole graph. In such a way all point-like interactions, leading to a unitary time 
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Figure 1: (Color online) A star graph F with scattering matrix in the vertex and all edges 
of length L. We consider the entanglement entropy of the edge i (red) with respect to all the 
others. 



evolution of the system, are covered. The most general boundary conditions, implementing this 
natural physical requirement, are [HI |12] at the vertex 

M 

[A(I - U),, ^j{t, 0) - t{I + lJ),,id,^j){t, 0)] = , (2.2) 

i=i 

where U is an arbitrary M x M unitary matrix and A a real parameter with the dimension of 
mass. To fully specify the problem we also need to impose boundary conditions at the external 
ends of the edges. The most general ones are 

{d,ij,){t,L)=fiiiJi{t,L), (2.3) 



where /ij are again real parameters with the dimension of mass. Eq. (2.3) is the familiar Robin 

(mixed) boundary condition. Notice that Eq. (2.2) extends this condition to the vertex V of 
the graph. 

It has been established in |Hlll2] that the point-like interaction, induced by (2.2), generates 
the scattering matrix 

.(jA- [A(I-U)-MI + U)] 

[A(I - U) + MI + U)] • ^^-^^ 

Besides of unitarity 

E{kp{k)=I, (2.5) 

E){k) satisfies Hermitian analyticity 

§\k) = 8{-k) (2.6) 

as well. Notice also that 

§(A)=U, §(-A)=U-\ (2.7) 



showing that the unitary matrix U, entering the boundary conditions (2.2), is actually the 
scattering matrix at the scale A. 

The main difficulty in solving the Schrodinger equation (2.1) on the graph F is the mixing 
between the different edges codified in the boundary conditions (2.2) and (2.3). In order to 
simplify this problem, we impose that the boundary conditions at the ends of each arm are all 
the same, in such a way to restore (at the level of the Hamiltonian) permutation symmetry of 
the edges of the graph. It should be clear from the physical point of view that, being interested 
in the thermodynamic limit with L,N ^ oo, the boundary conditions at L must not affect the 
final result. Thus we assume from now on that 



(2.8) 



Under this condition, Eqs. (|2.2|) and (|2.3|) can be rewritten in equivalent forms without mixing. 

alizi 



Indeed, let us introduce the unitary matrix U diagonalizing U, namely 

TT 



diag (e-2*"i,e-2^"^ 







TT 



< «j < 

2 2 



Remarkably enough, U diagonalizes also §(A;) for any k: 

k + ir]i k + iri2 



Sd{k) = U^S{k)U = diag 



k + irjM 



k ~ ir]i' k — ii]2' ' k — ii]M 



where 

rji = Atan(a;j) . 

It is quite natural at this point to introduce the fields 

M 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



which obviously satisfy Eq. ( 2.1[ ). In terms of the boundary conditions ( 2.2[2.3 ) decouple. 



lj,ipi{t,L) 



(2.13) 
(2.14) 



defining a simple spectral problem on the tensor product "H = <S>i=i ^^[O, L], which is analyzed 
below. 

It is worth stressing that ipi{t, x) is a superposition of the values of the original field ipi{t, x) 
at the same distance x from the vertex, but on different edges of the junction. Being so delocal- 
ized, (pi{t, x) is unphysical and provides only a convenient basis for dealing with the boundary 
conditions. The physical observables and correlation functions will be always expressed in terms 
of the physical fields ipiit.x). 



The eigenfunctions of —d^, obeying ( |2.13[ ) and ( |2.14[ ) are 

(j)i{k,x) = I 



k>0, 



(2.15) 
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where q are some constants to be fixed below and k satisfy 

\k — irji J \k + ifi 

In order to determine tlie spectrum of k explicitly, we simplify the problem further by requiring 
scale invariance. 



(2.16) 



2.2 The scale invariant case 

Scale invariance of the boundary conditions ( |2.13[ ) and (2.14) implies the values 





0, 
oo 



0) , Neumann b.c. 

7r/2) , Dirichlet b.c. 



oo [ai 

In other words, the critical points are fixed by the (M + 1) vector (/i, ?7i, 772, ?7j\/) 



(2.17) 
whose 



components take the values (2.17). For an edge i one has the following possibilities: 



(a) /i = (Neumann condition at x = L): Eq. (2.16) gives 



2ikL 



k + irji 
k - ir]i 



which gives 



Vi 



(t)i{n,x) 



00 



cos 



sm 



[n - 1)tt- 



X 

Li 



n = 1,2, 





x' 











n = 1,2,...; 



(2.18) 



(2.19) 
(2.20) 



(b) /i = 00 (Dirichlet condition at x = L): Eq. (2.16) implies 

' k + irji 



2ikL 



k - ir]i 



which gives 



Vi 
Vi 



(pi{n,x) 





x' 











00 



COS 



2 / X 

sin Inn 
Ij \ Ij 



n = l,2,... 
n = l,2,... 



(2.21) 



(2.22) 
(2.23) 



Notice that any of the sets (2.18, 2.19, 2.22, 2.23) represent a complete ortho-normal system 
in L2[0,L]. 
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2.3 Scale invariant scattering matrices 



Observing that the eigenvalue of § is 1 for rji = and —1 for rji = oo one concludes that the 
most general scale-invariant scattering matrix, compatible with a unitary time evolution, is 
given by 

§ = W§dW^ §d = diag(±l,±l,---±l), (2.24) 

where W is a generic M x M unitary matrix. From the group-theoretical point of view, any 
critical S matrix is a point in the orbit of some §d under the adjoint action of the unitary group 
U{M). Obviously, one can enumerate the edges in such a way that the first p eigenvalues of 
§ are and the remaining M — p are —1. The cases p = M and p = correspond to S = I 
and S = — I and are not interesting for the entanglement. In these two cases in fact, the single 
wires are decoupled (there is no transmission), which implies a vanishing entanglement. 



It follows from (2.24) that besides being unitary, § is also Hermitian (in agreement with 
(2.6) at criticality). Therefore, all diagonal elements Su are real. Notice however that in general 
S is not symmetric. If this is the case, time reversal invariance is broken [19]. 

In the nontrivial case p = 1 the most general 2x2 scale-invariant scattering matrix depends 
on two parameters and can be written in the form 

^) = ( 2ee-^' l -e' j ' ^ ^ ^' ^ ^ [0, 27r) . (2.25) 

Time reversal invariance is broken for 9 ^ 0, tt. 

For M = 3 one has two families corresponding to p = 1 and p = 2. In order to avoid 
cumbersome formulae, we display only two representatives of these families, namely 



=i(ei,e2) = , , , , A 262 -el + el-l 2e^e2 , (2.26) 



2ei 


262 




1 


- (-2 _ p2 
^1 ^2 


262 




1 




26i62 




26i62 








^1 ^2 ^ 


26162 






2ei 


26i62 


+ f2 _ 
^1 ^ ^2 


1 




262 


261 


262 




1 


- 6^ - f 2 
^1 ^2 



Sp=2(ei, 62) = ,^2^2 I 26162 "6? + 6^ - 1 262 | , (2.27) 

i ~r 6]^ + 6r 

where 61 2 G M. 



2.4 Two-point correlation function 

Now we are in position to construct the physical field ipi{t, x) and the relative two-point function 
needed in the computation of the entanglement entropy. First of all, we write the unphysical 
field in terms of the eigenf unctions (f)i{n,x) 

00 

<f,{t, x) = J2 e"'"'^"^* Hn, x)a,{n) , (2.28) 

n=l 
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where the fermion annihilation and creation operators satisfy standard anti-commutation rela- 
tions 

[ai{m) , a](n)]+ = SijSmn , [ai(m) , aj{n)]+ = [aj(m) , aj(n)]+ = , (2.29) 
and the energies are given by 

Uhz]^ i^p<^<M, 

for /i = and /x = oo respectively. Notice that different "unphysical" edges may have different 
dispersion relation, which is not a problem because these edges are totally isolated from each 
other. By means of (2.12) one gets the physical fields 

M M oo 

Ut,x) = J2^l^,{t,x) = 5^5^4e-'-^(")*0,(n,x)a,(n). (2.32) 

j=l j=l n=l 

One easily verifies that 

[^,(t,x) , ^J(t,y)]+ = 5i,5{x - y) , (2.33) 

which fixes the normalization of the fields. 

The equal time two-point correlation function of the physical field ipi{t,x) on a given state 
l^') is 



Cr^{x,y) ^ ml{t,x)ij,{t,y)\^) 

M oo 
k,l=l n,m=l 

where the correlator {'^\a]^{n)ai{m)\'^) can be deduced from the action of the algebra generated 
by {ai{m) , aj(n)} on the state In particular, we are interested in the case when is the 
ground-state of the system formed by fermions in the whole junction. It is then useful to 
rewrite N as 

= M^^ (2.35) 

where A/" represents the average number of particles for each wire. The action of the annihilation 
and creation operators on the ground-state is obvious since it is annihilated b y all aiim) with 
m> J\f and so {^\a\{n)ai{m)\^) = 5ki5nmQ{.-^ — n). Using this relation, Eq. (2.34), restricted 
to the same edge {i = j) which is needed actually for computing the entanglement entropy, 
becomes 

M Af 

C^{x, y) = J2Yl l^'^' I' Mn, x)Mn, y) , (2.36) 

k=l n=l 



where Af can also be interpreted as an ultraviolet cut-off for the series in (2.34). 
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It is convenient for what follows to rewrite (2.36) in more explicit terms. For this purpose 
we consider any critical point characterized by the integer 1 < p < M, i.e. a scale invariant 
scattering matrix with p eigenvalues equal to +1. The two sums in (2.36) factorize and one 
gets 



k=l 



n=l 



where 



f-{n,x) 

Because of (2.24), the sums over k in (2.37) give 

Em 





M N 
k=p+l n=l 




/ "j- cos \{n — 




-- 


= 0, 


/ 1; COS \{n — 




-- 


= oo , 


'l sin [(^- - 




^^ 


= 0, 


'l sin h^f] 
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= oo . 



(2.37) 



(2.38) 



(2.39) 



k=l 



M 

fc=p+i 



;i-§n) . 



(2.40) 



Moreover, using that S is both unitary and Hermitian, one has 

M 



(2.41) 



j=i 



where is the total transmission probability from the edge i to the rest of the graph. 

In conclusion, the correlator (2.36) can be fully expressed in terms of the transmission 
probability and the one-particle wave functions as follows 



C!Iix,y) 



1 - 



— \ ^ 

^ n=l 



{x,n)U{y,n) 



^ n=l 



{x,n)f^{y,n) 



(2.42) 

which is the basic input for deriving the entanglement entropy in the next section. We observe 
that involves only the diagonal elements of S and consequently, does not depend on the be- 
havior of S under transposition. Therefore, contrary to the conductance [H], the entanglement 
entropy in our case is not sensitive to the breaking of time-reversal invariance. 



3 Entanglement entropy 



In order to compute the bipartite Renyi entanglement entropies defined as in Eq. (1.1) of a 



subsystem A in the ground-state our star graph, we use the method recently introduced in 
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Refs. [381 [39] . The starting point to deal with a system made of a finite number of particles in 
continuous space is the Fredholm determinant 



i(A) = det [X6a{x, y) - Ca{x, ; 



(3.43) 



where CA{x,y) is the restriction of the correlation matrix C{x,y) defined in Eq. (2.34) to A, 
i.e. Ca = PaCPa, where Pa is the projector on A. The same definition holds for 6A{x,y) = 
PaS{x — y)PA- Following the ideas for the lattice model |43j, D^(A) can be introduced in such 
a way that it is a polynomial in A having as zeros the eigenvalues of Ca- Since we are dealing 
only with free fermions in the bulk, the reduced density matrix pA is Gaussian and so one 
can easily derive |13], |38l [39] 



In Trp° 



dX rflnDA(A) 



1 — a 

where the integration contour encircles the segment [0, 1], and 

1 



ea(A) 



1 — a 



In [A" + (1 - A)' 



(3.44) 



(3.45) 



For a — )■ 1, ei(A) = — xlnx — (1 — x) ln(l — x) and Eq. (3.44) gives the von Neumann entropy. 

The Fredholm determinant is turned into a standard one by introducing the reduced overlap 
matrix A (also considered in Ref. [15]) with elements 

dz(t)^{z)(t)^{z), n,m = 1,...,D, (3.46) 

Xl 



where in general (f)n{x) represent the eigenfuctions corresponding to the D lowest energy level 
which are occupied in the ground-state of the system with D degrees of freedom. The matrix 
A satisfies TrC^ = Tr A'^ and so [39] 



\nBA{X) 



gTrC^ 

k=l 



kX^ 



k=l 



TrA*^ 

Ta^ 



D 



In det [AI 



5^1n(A 



(3.47) 



m=l 



where am are the eigenvalues of A and D is its dimension to be specified later. Inserting (3.47) 



in the integral (3.44), we obtain 

D 

2m 



dX ea{X) 



m=l 



X 



D 
m=l 



a 



-Tr In [A" + (I 



(3.48) 



as a consequence of the residue theorem. 

In the following we will be interested only in the entanglement entropy of any edge i of the 
wire with respect to the rest of the junction in the global ground-state of the star graph. As we 
have seen above in Eq. (2.42), the two-point correlation function for finite number of particles 
in the full star graph, can be written in the form (we omit the edge index i hereafter) 

C''{x,y) = J2x{x,n)x{y,n), (3.49) 



n=l 
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where xi^^''^) ^ire proportional to the one-particle eigenfunctions f±{x,n) in Eq. (2.42). Then 
the Renyi entanglement entropy of the subsystem represented by a single wire is given by Eq. 
(3.48) where the eigenvalues of A are numerically calculated from the overlap matrix built 
with the correlation function above. 



3.1 The case M = 2 



It is instructive to consider first the case with two edges only. In this case, one has actually a 
segment of length 2L with a point-like conformal defect placed in the middle. This situation 
has been investigated on the lattice [351 ES] by other methods and represents a useful check for 
our framework. We set L = 1 for simplicity and consider the case /i = oo (Dirichlet condition 
at X = Combining the explicit form (2.25) of the §-matrix with (2.41), one obtains the 
following transmission amplitudess 



Ti = T2 



2e 

= T 



which are the square root of the transmission probability. Accordingly, 



2Af 



n=l 



with 



X{k,x) 



71^5^2 cos {kTTl) 

7Ti?V2sin(fc7rf) 



k = 1,3,. ..,27V- 1, 
k = 2,4,...,2Ar, 



(3.46), reads in our case 



(3.50) 



(3.51) 



(3.52) 



where N = 2M is the total particle number according to eq. (2.35). The matrix A, defined by 



A. 
A. 



1 



I 

-5„ 



:i + e2)""^"' 

2e 2n 
[1 + e^) 7r(ra^ — th?) 
2e 2m 



(1 + e^) 7r(m^ — n^) 



2\ ' 



m,n — odd - 



m — odd , n — even , 
m — even , n — odd . 



(3.53) 
(3.54) 
(3.55) 
(3.56) 



Since C 
by 



C22, the entanglement entropy of the edge 1 equals that of the edge 2 and is given 



N 



S^''\T;N) 



1 — a 



Tr In [A° + (I 



y^ea(an) 



(3.57) 



n=l 



"'^The case /i = (Neumann boundary conditions aX x — L) can be treated along the same lines. 
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an being the eigenvalues of the matrix (3.53 3.56). Using (3.57), we will show below that 



Si^){T-N) =C^''\T) lnA^ + 0(l) 



(3.58) 



with a pre-factor C^"\T) that depends on the transmission amplitude and not only on the 
central charge c = 1. 

Before considering a generic value of a, it is instructive to discuss the Renyi entropy a = 2. 



In this case Eq. (3.57) takes the simple form 



k=l 



Si^)(T- N) = -Tr In [I - 2E(T)] = V — TrE'=(T) 
where the combination 



(3.59) 



E(T)=A(I-A) (3.60) 
represents the natural variable for performing the computation. The main idea at this point 



is to reduce the evaluation of (3.59) to the case of full transmission T = 1, when the defect is 
absent and one can use therefore the known [10] behavior of Tr for free fermion gas on the 
interval, namely 

Tr E^T = 1) = — ^1^^ _ [J In iV + 0(1) . (3.61) 
For this purpose we first establish the fundamental relation 



TrE*=(T) = T^^ TrE''(r = 1) 



(3.62) 



which captures the impact of the point-like interaction in the junction on the entanglement 



entropy (3.59). In order to prove (3.62), we exploit the property that a reordering of rows and 
columns of E does not change the trace of E'^, we are interested in. We then write the overlap 
matrix in the following block form 



T] 



l+e- 



■I 



(3.63) 



with quadratic blocks of size M . Here I is the A/" x A/" identity matrix and with respect to 
(3.53 3.56) we have reordered the lines and rows of A in such a way that the upper block on 



the left and the lower one on the right coincide with (3.53) and (3.54) respectively (i.e. after 
this reordering the first lines and rows are the odd indices n, m, while the right-lower block 
is formed by even n,m). Bi and B2 are T-independent matrices, whose explicit form is not 



essential for the proof. Using the representation (3.63) and the relation (3.50), one gets 



E 




-T] 



T 



Ui 10)2 

I 



(3.64) 
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which proves ( [3^ . Finally, plugging (|3]6l]|3]62|) in one obtains 



1)!]^ 

^ 7T^k{2k-iy. 



T 



2k 



arcsm 



where we used 



2 arcsin x . 



We stress that each integer value of a > 2 can be treated in analogous way. For a 
finds for instance 



Ci^)(T'\ = \^ 3^[(A: - 1)!]^ 2k 



arcsin 



TVs 



(3.65) 



(3.66) 



3, 4 one 



(3.67) 



and 



k=l 



6n'^k{2k - 1)! 



2 

3^ 



arcsm^ ( ^y^ + v^i j + arcsin^ f^\/2 - ^2 



(3.68) 



3.1.1 Result for generic integer a 

In order to obtain the result for generic integer a we first need to write the combination 



A" + (I — A)"" in terms of the matrix E. This can be achieved by formally inverting Eq. (3.60) 
as 



'1 ± Vl - 4E) 



(3.69) 



with an ambiguity in the choice of the sign refiecting the degeneration of the spectrum of E. 
However, in the needed combination 



A" + (I - A)" = 2"" (1 ± Vl - 4E)° + {lT Vl - 4E)" 



(3.70) 



the choice of this sign is unimportant. Notice that although the apparent presence of a square- 
root, for integer a the above expression is a polynomial in E of degree [a/2j , that is the integer 
part of a/2. Using the binomial theorem (1 + x)" = X]fc=o (fc)"^'^' have 

A" + (I-A)° = 2-"^ r^Vl -4E)'=/2(i _ ^3 71) 

k=o ^ 

L"/2J . s 
fc=0 ^ ^ 



12 



which makes the polynomial form explicit. Using again the binomial theorem for (1 — 4E)^, we 
have 



L"/2J 



L«/2J 



fe=o V^'^/ p=o p=o 

where we defined 



L«/2J 

fe=p 



2A;) (p 



Q\ r(a-p)r(p + i/2) 

2pJ V^r(Q;) 



(3.73) 



In order to calculate the Renyi entropies, we expand in series of E the quantity 
In [A" + (1 - A)"], obtaining 



L"/2J 



L"/2J 



ln[A" + (1 - A)"] = ln[l - ^ Vp{-A¥.Y] = - XI " [ ^p^"^" 

p=l 3=1 p=l 



E) 



(3.74) 



Using now the multinomial identity we have 



L«/2J 



ln[A" + (l-A)1 = -5^-,5] 



n K(-4E)^]^'' (3.75) 



where we introduced the symbol Yl'k ^'^^ constrained sum ^^.^ with '^ki = j. We 



now introduce a sum over which will be equal to Y^pkp with the help of a Kronecker delta: 

L«/2J 



in[A" + (1 - A)«] = - E E E ^T^V;-T(-4E)"5i.,E..^. H ■ 

K=l j=l ki La/2J- p=i 



(3.76) 



Using the contour integral representation of the Kronecker delta over the unitary circle |2;| = 1 



1 

27ri 



dzz 



a-b-l 



(3.77) 



we have that the above expression equals 



1 

27H 



L«/2J 



n( 



1 

27ri 



oo L"/2J 

E(-4E)^ ^ (i^^^-Mn(l+ E V^) 

K=l p=l 



-p\kp 



K 

oo 



;i + Vi + ^T + (1 - ^1+^)" 



2" 



(3.78) 
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where, in the last hne, we recognized the expansion of the function from which we started from 
for a complex argument (no matrices). 

The integral over z can be performed with standard techniques of integrals on the complex 
plane. The various contributions come from the discontinuities at the cuts which are placed 
between the zeros of the argument of the logarithm, i.e. at the satisfying 



1 + + 



All the solutions of this equation are simply found as 

2^(2p-l) 



1 + 1 



0, 



cos 



2a 



with p = 1, . . . , a/2 . 



Thus the integral is given by 



1 
27ri 



dzz 



K-l 



In 



;i + ^TT^)" + (1 - 

2" 



L"/2J 



implying 



ln[A° + (1 - A)^ 



K 



-1 



K=l 



K 



L«/2J 

p=i 



cos 



2K 



TC- 



2p-l 
2a 



K 



(3.79) 



(3.80) 



(3.81) 



(3.82) 



In this final form, it is straightforward to take the trace using Eq. (3.62) to obtain 

2 L"/2J 



Trln[A" + (1 - A)"] = - At IniV V t^k^k HK 1)!] ^ 

K=l ^ ' p=l 



COS 



2K 



TT- 



2p 



2a 



(3.83) 



Inverting the order of the sums, the sum over K can be now performed using Eq. (3.66) and 
we have 

o L"/2J 



Trln[A" + (1 - A)"] = -InA^— ^ 
that leads to the coefficient 

C(")(T) 



arcsm 



Tcos 



(2p- l)7r 



2a 



a 



1_2}-^ 

p=l 



arcsm 



Tcos 



(2p - l)7r 
2a 



(3.84) 



(3.85) 



For a = 2,3,4 it coincides with the result reported in the previous subsection. 

We observe that C^"\0) = and C^°'\l) provide useful checks. The value T = describes a 
totally reflecting defect; the two edges are isolated and indeed the entanglement vanishes. For 
T = 1 we have 

L"/2J 



p=i 



IT (2p — l)7r 
2 2^ 



1 
12 



1 + 



a 



(3.86) 



This corresponds to full transmission, i.e. the defect is absent and one is considering the 
entanglement entropy of half of a system of length 2L with boundaries obtained in (3^ and 
compatible with the standard CFT result \27\ . 
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3.1.2 The analytic continuation 



In order to investigate non-integer values of a, we exploit the results of Ref. [351 [36] for the Ising 
and XX spin-chains. Using methods based on corner transfer matrix and conformal mappings, 
Eisler and Peschel derived the entanglement entropy of the subsystem on (let us say) the right 
of the defect, for a chain of length 2L with a defect in the middle. For the XX chain, that 
after a Jordan- Wigner transformation corresponds to a lattice gas of free spinless fermions, the 
result can be written in the form [35] 



5(")(r;L) 



C^l\t) lnL + 0(l) 



with C^""* (T) given by 



ct\T) 



a) 



dx In 



1 _^ Q~2aui{x,T) 



where 



uj{x, T) = acosh 



cosh(x) 
7^ 



(3.87) 



(3. 



(3.89) 



Actually only the result for a = 1 has been reported in Ref. [35] , but the derivation for general 
a is straightforward from the results reported there. A result of a = 2 has been reported in 
[SB] and for general integer a in [3B] . For a simple comparison between our work and Ref. [3B] , 
we mention that the transmission amplitude T in Ref. |i35j is called s. 

For integer a, it is straightforward to check numerically that C^^\t) and C'-°^(T) in Eq. 



(3.85) are equal, as it is possible to show analytically [46J. This coincidence does not come 
unexpected. Indeed, since Eq. (3.87) is valid for finite density N/2L = 1/2 on the lattice, it 
can be turned in the entanglement entropy as function of A^, simply by replacing L with A^, 
but assuming (as we did in on the basis of numerical data) that the dependence on T is 
universal. Taking now the continuum limit, we straightforwardly deduce that 



C(°)(T) =ci"^(T) 



(3.90) 



However, the above computation proves the universal T dependence with no assumption. The 
computation of Ref. [16] also shows that the result in Eq. (3.88) is the analytic continuation 
of Eq. (3.85) to non-integer values of a. In particular, in the von Neumann case {a = 1), the 



integral can be performed and one has 



Ci\T) 



{[(1 + T) ln(l + T) + (l-T) ln(l 



T)] InT + (1 + T)Li2(-T) + (1 - T)Li2(T)} . 

(3.91) 



3.1.3 Comparison with numerical computation 

We now turn to briefly present the numerical data for the entanglement entropies which have 
been fundamental for the conceptual understanding that led to the exact computation in the 
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Figure 2: (Color online) The von Neumann edge entanglement entropy in the two-wire junction 
for T = 3/5 and T = 4/5 where we considered only even values of A^. Left: We plot S^^\T, N) 
vs A^. Right: The subtracted quantity S^^\T,N) — C^^\T)\nN, where C^^\T) is given by 
Eq. (3.88). The lines show fits of the data for > 400, where the correction to the leading 
behavior C^^\T)\nN is a polynomial bo + bi/N + 62/^^- 



previous subsection. Furthermore numerical computations give non trivial insights about the 
corrections to the asymptotic behavior. In particular we anticipate that in analogy to what 
found for the half-space Renyi entanglement entropies in homogeneous systems [39l W7\ . the 
leading suppressed corrections turns out to be 0{N~^^°'). 

The numerical estimates of the factor C^^\T), obtained from our results for the entanglement 
entropy up to ~ 500, perfectly match the function (3.91), within a precision better than 
O{10-^). Fig. |2| shows the results for T = 3/5 and T = 4/5 for even A^. Fits to the T = 4/5 
data for 400 < A^ < 500 



S^^^ (T; A^) = a In A^ + 60 + fei/A^ + b2/N^ + b^/N^ 



(3.92) 



give a = 0.118841, 69 = 0.342056, 61 = —0.1404, etc., where a should be compared with 
C(i)(4/5) = 0.11884065.... Analogous resuhs are obtained for T = 3/5; we find a = 0.076078 to 
be compared with C^^\3/5) = 0.07607750. 

A complete agreement is also found for the Renyi entropies. For a = 2 for instance, the 
data fit the asymptotic behavior 



S^'^\T; N) = C^'^\T) In A^ + 60 + bi/N^''^ + 62 /A^ + b^/N^^'^ + 



(3.93) 



where the leading coeffi cient C^^^T) is given by Eq. (3.65), as shown by Fig. [3] 

In both Eqs. (3.92) and (3.93) we use for the subleading corrections the form 0{N~^/°') 



We mention that considering also odd values of the particles number A^, we observe also in 
the graph entanglement entropies corrections to the scaling which depends on the parity of A^, 
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Figure 3: (Color online) The a = 2 Renyi edge entanglement entropy S^'^\T, N) in the two-wire 
junction for T = 3/5 and T = 4/5 for even number of particles N. The line show the curve 
C(2)(T) \nN + bo + bi/N'^^^ + b2/N where the coefficients hi are fitted using the data for > 400. 



in analogy to what observed in the absence of defects both in the continuum [SH] and on the 
lattice H?]. 



3.2 The case of M > 2 edges 



The generalization to M > 2 edges is now straightforward. The novelty is that now the 
entanglement entropy on the edge i with respect to the whole junction depends on i. As 



expected from (2.42), this dependence is encoded in the transmission amplitudes Tj, which 



for M > 2 do not coincide in general for different i. In order to investigate this aspect, it is 
convenient to introduce the variable 



^-2 



l-T; 



(3.94) 



The two point-function (2.42) takes now the form 



n=l 



n=l 
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Equivalently, using (2.38|2.39), one finds in the case /i = oo 



2N 



(3.96) 



n=l 



with 



'Tiv^cos (A;7rf) 



k = 1,3,...,2A/' - 1, 
k = 2,4,...,2Ar. 



T^Tfx/2sin(fc7rf) , 

The reduced overlap matrix A now also carries an edge index i = 1 , . . . , M and is given by 



(3.97) 



mn 



2(1 - r'.)s„ 



mn 



mn 



2T, a/1 - T? 



2n 



2T,a/1 - T: 



7r(r?,^ — m^) 
2m 

7r(m^ — n^) ' 



m,n — odd , 
m,n — even , 

m — odd , n — even . 
m — even , — odd , 



(3.98) 
(3.99) 

(3.100) 
(3.101) 



Comparing the matrices (3.53 3.56) and (3.98 3.101) one concludes that the entanglement en- 
tropy S'f"'' of the edge i with respect to the whole junction depends on the edge via the trans- 
mission amplitude 

T, = 2T, ^l-T2^y/l-§2 (3.102) 

and can be expressed by 

^J"^(Ti; N) = S^''\T,; N) , (3.103) 



where S^°'^ is the entropy (3.57) in the case M = 2. Therefore, the asymptotic behavior for 
large is given by 

Si"\Ti;N) =C(°)(Ti)lnA^ + 0(l), (3.104) 



where C*-"-* is the universal function defined by (3.85|3.88, 3.89). At this point all the numerical 



results in the previous subsection concerning C^"^ apply. 



4 Schrodinger junction with harmonic potential 

In this section we consider a star graph F with M infinite edges and a harmonic potential 
V{x) = \moj'^x'^ trapping the gas in the bulk (harmonic trap). The Schrodinger field ipiit^x) 
thus satisfies 

idt + —dl--mu'^xAiJi{t,x) = x > . (4.105) 

2m 2 / 

Since V[x) defines a self-adjoint multiplication operator, the vertex boundary conditions con- 



trolling the self-adjointness of the total Hamiltonian —d^ + V{x) are still given by (2.2). Ac 



cordingly, the critical S matrices are parametrized by (2.24). The field (pi{t,x), defined by 
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(2.12), satisfies 



(5,¥P,)(t,0) = 0, l<t<p, 
(^.(t,0) = 0, p<i<M, 

for all t. The eigenfunctions of —d^ + V{x), obeying these boundary conditions, are 



f+{n,x), l<i<p, 
f-{n,x) , p < i < M , 

where f±{n, x) are expressed in terms of the Hermite polynomials as follow^ 

1 „ , . _.2, 



f-{n,x) 



7ri/4^22"-i(2n)! 
1 



i^2n(x)e-^ 

x)e-^'/2 



2n+H 



n = 0, 1,, 



0,1, 



7ri/4^22"(2n + 1)! 

Inserting ( |4l09pll0| in ( [3^ , one gets the two-point correlation function 



(4.106) 
(4.107) 



(4.108) 

(4.109) 
(4.110) 



2Af 



(4.111) 



n=0 



with 



=W,Hk{x)e- 




A; = 1,3, 



2A^. 



(4.112) 



= 0,2,...,2Af. 



At this point one can derive the relative A-matrix and compute the entanglement entropy. 
The numerical data, displayed in fig. 4 , confirm that in the harmonic case SriTi; N) has pre- 



cisely the behavior described by equation (3.104). An educated guess for the large- asymptotic 



behavior of the edge entanglement entropy is that it is the same as that of the hard-wall case, 
i.e. 

5(")(T; A^) = C("HT) lnN + 0(1), (4.113) 



where the functions C^"-* are the same as those in Eq. (3.88). Difference occurs at the level of 



the 0(1) term, as in the case without defect |3Hl HHl US]. 

In order to check it for nontrivial defects, we have computed the entanglement entropy up 
to A^ ~ 50 for the two-edge problem, and for two values of T, T = 3/5 and T = 4/5. The 
results for the von Neumann entropy are shown in Fig. |4j They are perfectly consistent with 
the conjecture (4.113). For example a fit of the data for T = 4/5 and even A^ to 



5^^) (T; A^) = a In A^ + 6o + fei/A^ + bi/N^ + h/N^ 



(4.114) 



give a = 0.118842, bo = 0.424421, etc., where a should be compared with C^^\T = 4/5) 
0.11884065.... 



^We set for simplicity ui = 1 and m — 1. 



19 



0.25 



1 ' 1 


1 

T=4/5 


- S'^'(T;N) - d^\T)lnN 


- 

- 




T=3/5 


1 


harmonic trap 

1 



0.00 0.02 0.04 0.06 0.08 



1/N 



Figure 4: (Color online) The von Neumann edge entanglement entropy in the two- wire junction 
with an external harmonic potential for T = 3/5 and T = 4/5. We plot the subtracted quantity 
S^^\T; N) - 6^\T) In A^, where C^^\T) is given by Eq. ( |3.88[ ). The lines show fits of the data 
for > 400 to a polynomial 6o + &i/iV + &2/A^^- 



5 Outlook and conclusions 



We have studied the entanglement entropies of one edge % with respect to the rest of a junction 
with M edges. Our main result is that when working with a finite number of particles iV in 
edges of finite length L, the Renyi entanglement entropies can be derived analytically, obtaining 

5^") =C(")(r)lnAr + 0(ArO), (5.115) 

where T is the total transmission probability from the edge i to the rest of the graph and the 
pre-factor is independent on the number of edges. The analytical computation of C*-"^ is given 



by Eq. (3.85) for integer a and its analytical continuation to non-intger a is given by the 



prefactor of Eisler and Peschel [35] obtained for the spatial entanglement of a line with a defect 



reported in Eq. (3.88). Clearly the value of the total transmission does depend on the number 
of edges and of the kind of junction. The same asymptotic behavior in also describe the 
entanglement entropies of systems in which on each arm there is a confining parabolic potential. 

We can turn the above asymptotic behavior in a more standard expression for the depen- 
dence of the entanglement entropies on the length of the subsystem (i.e. the edge in our case). 
Indeed, assuming a uniform density of particles we have N oc L and so 

S(") = C^") (T) In L + 0{L°), (5.116) 
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that is expected to be valid also for lattice models and in particular for M XX spin chains of 
length L joined at a single common vertex. 



Finally, we want to mention that the asymptotic result (5.116) has not yet been derived 
from conformal field theory. Although it is clear that such derivation must be possible because 
the CFT encodes all the required ingredients, the practical calculation is very cumbersome, as 
the similar (but different) result of Ref. [33] shows. 

Acknowledgements. 

We are very grateful to Ingo Peschel for an intense correspondence about an independent and 



different way to obtain Eq. (3.85) and for sharing the results of Ref. [46j before publication. 
PC research was supported by ERC under the Starting Grant n. 279391 EDEQS. He also 
acknowledges the kind hospitality at the Institute Henri Poincare, Paris. 

References 

[1] C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 1220 (1992); Phys. Rev. B 46, 15233 
(1992). 

[2] 1. Safi and H. J. Schulz, Phys. Rev. B 52, R17040 (1995). 

[3] C. Nayak, M. P. A. Fisher, A. W. W. Ludwig, and H. H. Lin, Phys. Rev. B 59, 15694 
(1999). 

[4] 1. Safi, P. Devillard, and T. Martin, Phys. Rev. Lett. 86, 4628 (2001). 

[5] J.E. Moore and X.-G. Wen, Phys. Rev. B 66, 115305 (2002). 

[6] H. Yi, Phys. Rev. B 65, 195101 (2002). 

[7] S. Lai, S. Rao, and D. Sen, Phys. Rev. B 66, 165327 (2002). 

[8] C. Bachas, J. de Boer, R. Dijkgraaf, and H. Ooguri, JHEP 0206, 027 (2002). 

[9] S. Rao and D. Sen, Phys. Rev. B 70, 195115 (2004). 

[10] C. Chamon, M. Oshikawa, and I. Affleck, Phys. Rev. Lett. 91, 206403 (2003); M. Oshikawa, 
C. Chamon, and I. Affleck, J. Stat. Mech. P02008 (2006). 

[11] T. Enss, V. Meden, S. Andergassen, X. Barnabe-Theriault, W. Metzner, K. Schonhammer, 
Phys. Rev. B 71, 155401 (2005); X. Barnabe-Theriauh, A Sedeki, V. Meden, K. Schon- 
hammer, Phys. Rev. Lett. 94, 136405 (2005); X. Barnabe-Theriauh, A. Sedeki, V. Meden, 
K. Schonhammer, Phys. Rev. B 71, 205327 (2005). 



[12] D. Friedan, cond-mat/0505084 cond-mat/0505085 



[13] S. Das, S. Rao, D. Sen, Phys. Rev. B 74 (2006) 045322. 



21 



[14] B. Bellazzini and M. Mintchev, J. Phys. A 39 (2006) 11101. 

[15] B. Bellazzini, M. Mintchev and P. Sorba, J. Phys. A 40 (2007) 2485. 

[16] C.-Y. Hou and C. Chamon, Phys. Rev. B 77 (2008) 155422. 

[17] S. Das and S. Rao, Phys. Rev. B 78 (2008) 205421. 

[18] B. Bellazzini, P. Calabrese and M. Mintchev, Phys. Rev. B 79 (2009) 085122 
[19] B. Bellazzini, M. Mintchev and P. Sorba, Phys. Rev. B 80 (2009) 245441, 
[20] A. Soori and D. Sen, Europhys. Lett. 93 (2011) 57007. 

[21] L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Rev. Mod. Phys. 80 (2008) 517; J. Eisert, 
M. Cramer, and M. B. Plenio, Rev. Mod. Phys. 82 (2010) 277; P. Calabrese, J. Cardy, 
and B. Doyon Eds, J. Phys. A 42 (2009) 500301. 

[22] J. Eisert and M. Cramer, Phys. Rev. A 72 (2005) 42112; I. Peschel and J. Zhao, J. Stat. 
Mech. (2005) PI 1002; R. Orus, J.I. Latorre, J. Eisert, and M. Cramer, Phys. Rev. A 73 
(2006) 060303. 

[23] P. Calabrese and A. Lefevre, Phys. Rev. A 78, 032329 (2008). 
[24] J. Cardy, J. Stat. Mech. (2010) P10004. 

[25] C. Holzhey, F. Larsen, and F. Wilczek, Nucl. Phys. B 424 (1994) 443. 

[26] G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, Phys. Rev. Lett. 90 (2003) 227902. J. I. 
Latorre, E. Rico, and G. Vidal, Quant. Inf. Comp. 4 (2004) 048 . 

[27] P. Calabrese and J. Cardy, J. Stat. Mech. (2004) P06002. 

[28] P. Calabrese and J. Cardy, J. Phys. A 42 (2009) 504005. 

[29] G. Levine, Phys. Rev. Lett. 93 (2004) 266402. 

[30] I. Peschel, J. Phys. A: Math. Gen. 38 (2005) 4327. 

[31] J. Zhao, I. Peschel, and X. Wang, Phys. Rev. B 73, 024417 (2006); 
J. Ren, S. Zhu, and X. Hao, J. Phys. B 42 (2009) 015504. 

[32] F. Igloi, Z. Szatmari, and Y.-C. Lin, Phys. Rev. B 80 (2009) 024405. 

[33] K. Sakai and Y. Satoh, JHEP 0812: 001 (2008). 

[34] E. S. Sorensen, M.-S. Chang, N. Lafiorencie, and I. Affleck, J. Stat. Mech. (2007) LOlOOl; 
J. Stat. Mech. (2007) P08003; E. S. Sorensen, N. Lafiorencie, and I. Affleck, J. Phys. A 
42, 504009 (2009). 



22 



[35] V. Eisler and 1. Peschel, Ann. Phys. (Berlin) 522 (2010) 679. 
[36] I. Peschel, 1109.0159. 

[37] V. Eisler and S. Garmon. Phys. Rev. B 82, 174202 (2010). 

[38] P. Calabrese, M. Mintchev, and E. Vicari, Phys. Rev. Lett. 107 (2011) 020601. 

[39] P. Calabrese, M. Mintchev, and E. Vicari, J. Stat. Mech. (2011) P09028. 

[40] P. Calabrese, M. Mintchev, and E. Vicari, 1111.4836. 

[41] V. Kostrykin and R. Schrader, Fortschr. Phys. 48 (2000) 703. 

[42] M. Harmer, J. Phys. A 33 (2000) 9015. 

[43] B-Q Jin and V.E. Korepin, J. Stat. Phys. 116 (2004) 79. 

[44] I. Peschel and V. Eisler, J. Phys. A 42 (2009) 504003. 

[45] I. Klich, J. Phys. A 39, L85 (2006). 

[46] I. Peschel and V. Eisler, 1201.4104. 

[47] N Laflorencie, E S Sorensen, M-S Chang, and 1 Affleck, Phys. Rev. Lett. 96, 100603 (2006); 
P. Calabrese, M. Campostrini, F. Essler, and B. Nienhuis, Phys. Rev. Lett. 104, 095701 
(2010); P. Calabrese and F. H. L. Essler, J. Stat. Mech. (2010) P08029; M. Fagotti and P. 
Calabrese, J. Stat. Mech. P01017 (2011). 

[48] M. Campostrini and E. Vicari, J. Stat. Mech. (2010) P08020; J. Stat. Mech. (2011) E04001; 
Phys. Rev. A 81 (2010) 063614. 

[49] M. Campostrini and E. Vicari, Phys. Rev. A 82 (2010) 063636. 



23 



